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1 INTRODUCTION 

In this paper, we review the Cauchy problem for the space-time fractional 
partial differential equation, which is obtained from the standard diffusion 
equation by replacing the second-order space derivative with a Riesz-Feller 
derivative of order a £ (0, 2] and skewness 9 {\6\ < min{a, 2 — a}), and the 
first-order time derivative with a Caputo derivative of order (3 G (0, 2] . 



The fundamental solution (Green function) for the Cauchy problem is 
investigated with respect to its scaling and similarity properties, starting 
from its Fourier-Laplace representation. In the cases {0 < a < 2 , (3 = 
1} and {a = 2 , < (3 < 1} , the fundamental solution is known to be 
interpreted as a spatial probability density function evolving in time, so we 
talk of space-fractional diffusion and time-fractional diffusion, respectively. 
Then, by using the Mellin transform, we provide a general representation 
of the Green functions in terms of Mellin-Barnes integrals in the complex 
plane, which allows us to extend the probability interpretation to the ranges 
{0 < a < 2 , < 13 < 1} and {1 < (3 < a < 2}. 

Furthermore, from this representation it is possible to derive explicit 
formulae (convergent series and asymptotic expansions), which enable us 
to plot the spatial probability densities for different values of the relevant 
parameters a, 9, (3 . 

2 THE GREEN FUNCTION 

By replacing in the standard diffusion equation 

—u{x,t) = -^u(x,t), -oo<x<+oo, t>0, (1) 

where u = u(x, t) is the (real) field variable, the second-order space 
derivative and the first-order time derivative by suitable integro- differential 
operators, which can be interpreted as a space and time derivative of 
fractional order, we obtain a sort of "generalized diffusion" equation. Such 
equation may be referred to as the space-time-fractional diffusion equation 
when its fundamental solution (see below) can be interpreted as a probability 
density. We write 

t D^u{x,t) = x Dq u(x,t) , -oo<x<+cxd, t>0, (2) 

where the a , 9 , (3 are real parameters restricted as follows 

0<a<2, \e\ < min{a,2 - a}, < (3 < 2 . (3) 

In Eq. (2) x Dq is the space-fractional Riesz-Feller derivative of order a 
and skewness 9 , and tD* is the time-fractional Caputo derivative of order 
(3 . The definitions of these fractional derivatives are more easily understood 
if given in terms of Fourier transform and Laplace transform, respectively. 



2 



For the space-fractional Riesz-Feller derivative we have 

(4) 

^(«) = | K |« e <sign k)^^ 

where /(k) = J 7 {/(x); ^} = f-^ e' lKX f(x) dx . In other words the symbol 
of the pseudo-differential operatoJi] x Dg is required to be the logarithm 
of the characteristic function of the generic stable (in the Levy sense) 
probability density, according to the Feller parameterization [H [2] , see also 
Refs. 0|. For a = 2 (hence 9 = 0) we have 3o(«) = = H^) 2 , 
so we recover the standard second derivative. More generally for 9 = we 
have J)%{k) = -\n\ a = -{n 2 ) a l 2 



so 

/ d 2 \ a/2 

^0=" -v. ■ (5) 



dx 2 

In this case we call the LHS of Eq.(5) simply the Riesz fractional derivative 
operator of order a . Assuming q / 1,2 and taking 9 in its range, one can 
show that the explicit expression of the Riesz-Feller fractional derivative 
obtained from Eq. (4) is 

x D%f(x) :=- [c + (a,9) x D«+c„(a,9) x D a _] f(x), (6) 

where 

. „. sin [(a - 9) tt/2] sin [(a + 9) vr/2] 

sm (Q7TJ sm(a7rj 
and the X -D± are Weyl fractional derivatives defined as 

d 



dx 



*I± a f(x) 



if < a < 1 



e £>£/0«0H d2 (8) 



Jt a f(x) 



if 1< a < 2 . 



(ix 2 

In Eq. (8) the X I± (fi > 0) denote the Weyl fractional integrals defined as 

1 '" ""-'/(Ode, 

Ox > 0) (9) 




Let us recall that a generic linear pseudo-differential operator A, acting with 
respect to the variable x g JR , is defined through its Fourier representation, namely 
J+~e iKa: yl[/(a;)]da; = A(/s)/(/s), where is referred to as symbol of ^4 , given as 

= (Ae- iKX ) e +iKX . 
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In the particular case 6 = we get c+(a,0) = c_(a,0) = l/[2 cos (an/2)] , 
and, by passing to the limit for a — > 2~ , we get c+(2, 0) = c_(2, 0) = —1/2 . 
For a = 1 we have 



with X D f(x 



x Dlf(x)= [cos(#7r/2) x Z)i+sin(#7r/2),L>] f(x), 
d 



(10) 



ri.r 



/(x) , and 



-±[ x Hf(x)}, x Hf(x) = ± 



f(0 



(ii) 



In (11) the operator x i7 denotes the Hilbert transform and its singular 
integral is understood in the Cauchy principal value sense. 

The operator X D$ has been referred to as the Riesz-Feller fractional 
derivative since both Marcel Riesz and William Feller contributed to its 
definitional 

Let us now consider the time-fractional Caputo derivative. Following 
the original idea by Caputo [5] , see also [6j [8] , a proper time fractional 
derivative of order f3 S (m — l,m] with m G M , useful for physical 
applications, may be defined in terms of the following rule for the Laplace 
transform 

m— 1 

£{ t DPf(t);s} = sPf(s)-Y,^- 1 - k f {k) (0 + ), m-K(3<m. (12) 

k=0 

where f(s) = £{f(t);s} = / °° e — s * f(t) dt . Then the Caputo fractional 
derivative of f(t) turns out to be defined as 



tDSf(t) :- 



1 



/( m )(r)dr 



r(m - p) Jo (t - T y +l - m 

—f(t), (3 = m. 

dt m 



In other words the operator tD* is required to generalize the well-known 
rule for the Laplace transform of the n-th order derivative of a given (causal) 



m — 1 < f3 < m . 



(13) 



2 Originally, in the late 1940's, Riesz [9] introduced the pseudo-differential operator 
x Iq whose symbol is \K,\~ a , well defined for any positive a with the exclusion of odd 
integer numbers, afterwards named the Riesz potential. The Riesz fractional derivative 
x Dq := — xIq" defined by analytical continuation was generalized by Feller in his 1952 
genial paper [I] to include the skewness parameter of the strictly stable densities. 
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function keeping the standard initial value of the function itself and of its 
derivatives up to order n — 

The fundamental solution (or the Green function) G e a g(x,t) of Eq. (2) 
is the solution corresponding to the initial condition u(x,0 + ) = 5(x) . We 
note that if 1 < (5 < 2 we add the condition Ut(x, + ) = . 

In the particular case of the standard diffusion equation (1) the Green 
function is nothing but the Gaussian probability density function with 
2t , namely 



variance a 



G\ i(x, t) 



1 



20F 



-1/2 p -x 2 /(4t) 



(14) 



We note that in the limiting case {a = (3 = 2} we recover the D'Alembert 
wave equation with Green function G^ 2 { x i = [$( x + t)+ 6(x — t)]/2 . 

In the general case the application of the transforms of Fourier and 
Laplace in succession to Eq (2) yields 



-€^)G%{k,s) = S PG%{k,s)-sP-\ 



(15) 



3 The reader should observe that the Caputo fractional derivative differs from the usual 
Riemann-Liouville fractional derivative which, defined as the left inverse of the Riemann- 
Liouville fractional integral, is here denoted as tD 13 f(t) . We have, see e.g. [11] . 

I /- f(r)dr 



tITf(t) := 



r(m-j3) J (jt-r)P+i- 



m — 1 < (3 < m , 
P = m. 



When f3 is not integer and both fractional derivatives exist we have between them the 
following relation, see e.g. [7], 



tDSffl = tD p 



fc! 



1 < 8 < m, 



k-p 



m— 1 

^/(«) = ^/(*)-E/ (W ( fff )r(* -/j + i) 



1 < 8 < m. 



The Caputo fractional derivative, originally ignored in the mathematical treatises, 
represents a sort of regularization in the time origin for the Riemann-Liouville fractional 
derivative and satisfies the relevant property of being zero when applied to a constant. 
For more details on this fractional derivative we refer the interested reader to Gorenfio 
and Mainardi [7] and Podlubny [S]. 

Added Note: Nowadays the reader can find an exhaustive treatment of fractional 
derivatives in the treatise by A. A. Kilbas, H.M. Srivastava and J.J. Trujillo, Theory and 
Applications of Fractional Differential Equations, Elsevier, Amsterdam (2006). 
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namely 



G *>' S ' = ^T«W' (16) 

By using the known scaling rules for the Fourier and Laplace transforms, 
following the arguments by Mainardi et al. [TO] one can prove without 
inverting Eq. (16) that G a p(x,t) has the scaling property 

G%{x,t) = t~^K%(x/t^) . (17) 

Here xjt^l a acts as the similarity variable and K^g(-) as the reduced 
Green function. For the analytical and computational determination of the 
reduced Green function we can restrict the attention to x > because 
of the symmetry relation K^p{— x) = K~g(x). Extending the paper by 
Gorenflo et al. [12] where the space-time fractional diffusion equation has 
been considered with the restriction {1 < a < 2, 9 = 0, < (3 < 2}, 
Mainardi et al. [TO] have first inverted the Laplace transform and then have 
obtained the following Fourier integral representation of the reduced Green 
function 

•tfjn) dK. (18) 



I r+co 



K»Jx) = ^ I e~ iKX E, 



a > pK ' = 27 / " 



Here Ep denotes the entire transcendental function, known as the Mittag- 
Leffler function of order (3 , defined in the complex plane by the power series 

oo ri 

E ^^w^ry 3>0 ' <19) 

For detailed information on the Mittag-Leffler function the interested reader 
may consult e.g. [13J (Vol. 3, Ch. 18, pp. 206-227) and EH El . By 
using the convolution theorem for the Mellin transforrq3, the authors in Ref. 



4j f 

p+oo 



M{f(r);s} = f*(s)= /(rjf- 1 *-, 7 i<»(s)<72 



denotes the Mellin transform of f(r) , the inversion is provided by 

1 /"J+ioo 

X" 1 {.r( S );r} = /(r) = — / f(s)r~ s ds 

J y — ioo 

where r > , j = ?ft (s) , 71 < 7 < 72 • The Mellin convolution formula reads 

oo 

Kr) = J - p f(p)g(r/p)dpH h" (s) = f (s) g" (s) . 
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|10| have provided the Mellin-Barnes integral representation^ for the general 
case 0<a<2, 0</3<2asin Eqs. (2) and (3): 



ax2vri7 7 -ioo r(l - ^-s) Y(ps) T(l - ps) 2 a 

(20) 

where < 7 < min{a, 1}. Then, after changing s into — s and using the 
known property T(l + z) = , we can write the Mellin transform of 

f°V tos«cfc- r(i- a /a)r(i + a /q)r(i + 8 ) 

io )x dx - /, r(i- ps )r(i + p S )r(i + /3,/«)' ( } 

where — min{a, 1} < 5ft(s) < a . In particular we find f^ co K^Jx)dx = p 
(with p = 1/2 if = 0). We note that Eq. (21) is strictly valid as soon as 
cancellations in the "gamma fraction" at the RHS are not possible. Then 
this equation allows us to evaluate (in ]Rq) the (absolute) moments of order 
5 for the Green function if — min{a, 1} < 5 < a . In other words, it states 
that K e a p(x) = O ^x~( Q+1 )^ as x — ► +00 . When cancellations occur in 
the "gamma fraction" the range of 5 may change. An interesting case is 
{q = 2, 9 = 0, < (3 < 2} (time-fractional diffusion including standard 
diffusion), where Eq. (21) reduces to 



-OG 



1 ra + g'/q - !R<S)> - 1 - <22) 



This result proves the existence of all moments of order 5 > — 1 for the 
corresponding Green function which reads 

iifS fl (x) = -M g/2 (x) = / , 1 x s ds, x>0, (23) 



where Mp / 2 denotes a function of the Wright type introduced by Mainardi 
[T6| PT7] . see also [H [TBI IS]- F° r the case a = /3 (that we call neutral- 



The names refer to the two authors, who in the beginning of the past century 
developed the theory of these integrals using them for a complete integration of the 
hypergeometric differential equation. However, as pointed out in [13] (Vol. 1, Ch. 1, 
§1.19, p. 49), these integrals were first used by S. Pincherle in 1888. For a revisited 
analysis of the pioneering work of Pincherle (Professor of Mathematics at the University 
of Bologna from 1880 to 1928) we refer to the recent paper by Mainardi and Pagnini [15] . 
As a matter of fact this type of integrals turns out to be useful in inverting the Mellin 
transform. 
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fractional diffusion) we obtain from (20) an elementary expression 



ax2mJ 1 -ioo T(ps)T(l — ps) 
1 i p+ioo gin^ps) i x^sudf (a-0)] 



(24) 



ax 2ni J-y-ioo sin(-7rs/a) it 1 + 2x a cos[^(a — 9)] + x 2a 

where < 7 < a . 

We note that the Merlin-Barnes integral representation allows us to 
construct computationally the fundamentals solutions of Eq. (2) for any 
triplet {a,/3,6} (a 7^ 0) by matching their convergent and asymptotic 
expansions, as shown in [10J. Readers acquainted with Fox H functions 
can recognize in Eq. (20) the representation of a certain function of this 
class, see e.g. [U [UJ EO], [21] . Unfortunately, as far as we know, computing 
routines for this general class of special functions are not yet availabl^]. 

Let us now point out that, in the peculiar cases of space-fractional- 
diffusion, time- fractional- diffusion, and neutral- fractional- diffusion, non- 
negativity of the corresponding Green functions can be proven; being 
normalized in M these functions can indeed be interpreted as probability 
densities. In particular, for (3 = 1 and < a < 2 (strictly space-fractional 
diffusion) we recover the class of the strictly stable (non-Gaussian) densities 
L e a (x) exhibiting fat tails (with the algebraic decay oc |x|~( Q+1 )) and infinite 
variance, whereas for a = 2 and < (3 < 1 (strictly time-fractional 
diffusion) the class of the Wright-type densities Mp/ 2 ( x )/2 exhibiting 
stretched exponential tails and finite variance proportional to t@ . 

The meaning of probability density can be extended under proper 
conditions to the reduced Green function of the space-time-fractional partial 
differential equation (2) in virtue of the following identity, proven in Ref. 

□13, 



" / V _1 ^(n] L e a ( X /o ^, o<p<i, 

J o Mp/a (C)N e a (x/0j, 0</3/a<l. 

Then, due to the previous discussion, in the cases {0<a<2,0</3<l} 
and {1 < (3 < a < 2} we obtain a class of probability densities (symmetric or 



Added Note Nowadays the reader can find the representation of the fundamental 
solutions in terms of Fox H functions in the paper by F. Mainardi, G. Pagnini and 
R.K. Saxena, Fox H functions in fractional diffusion, J. Computational and Applied 
Mathematics 178, 321-331 (2005). 
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non-symmetric according to 9 = or 9 ^ 0) which exhibit fat tails with an 
algebraic decay oc . Thus, they belong to the domain of attraction 

of the Levy stable densities of index a and can be referred to as fractional 
stable densities. 



3 CONCLUSIONS 

The above analysis of the Cauchy problem for the space-time fractional 
diffusion equation shows the fundamental and combined roles of the integral 
transforms of Fourier, Laplace and Mellin type. By the aid of the transforms 
of Fourier and Laplace the scaling and similarity properties of the Green 
function can be easily derived. The Mellin transform allows us to obtain 
for the Green function a general representation in terms of Mellin-Barnes 
integrals (hence a computational form in terms of convergent and asymptotic 
series), and the extension of its probability interpretatiorjll 
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